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Abstract 

Discrete truncated power is very useful to study the number of nonnegative 
integer solutions of linear Diophantine equations. In this paper, by using the Laplace 
transform and the theory of toric arrangement, we show that discrete truncated 
power is a periodic piecewise polynomial on the shifted integral lattice cone. Based 
on the toric reduction method in the real field, we give a toric arrangement method 
to compute discrete truncated power. 
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1 Introduction 



Let be the collection of s- dimensional vectors whose components are 
integers, and R+ be the collection of all nonnegative real numbers. Let X be 
a multiset (i.e. the set whose elements can be same) of n vectors 
where Oj G Z''\{0}, i = 1, . . . ,n, such that for xi, . . . , x„ G R+, the equation 
J27=i ^i'^i = does not have nonzero solutions with respect to Oj, i = 1, . . . , n. 
Thus we see that 

n 

tx{x) := = (A, . . . , G N" : ^ f3,a, = x}, (1) 

1=1 

is finite for every x G Z**, where denotes the cardinality of the set A, 
and N" denotes the collection of n-dimensional vectors whose components are 
nonnegative integers. Then tx is called the discrete truncated power or the 
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partition function. Discrete truncated power was first introduced by Dahmen 
and Micclielli [9]. Tliey exploit the relationship between discrete truncated 
power and multivariate splines, and propose some important properties of 
discrete truncated power [lOj. In particular, they present a recursive method 
to compute discrete truncated power, which is based on the following recursive 
property: 

oo 

tx{a) = '^tx\{a,}{a - jai), i = l,...,n. 

j=0 

On the other hand, Laplace transform is found to be a useful tool to study 
the discrete truncated power. To see this, we let 

n 

Ilx{ci) '■= {x G R"| ^ Xitti = a, Xi> 0} 

i=l 

be a polytope. Then the value of tx{(y) is the number of integral points in 
Ilx{oi)i and the Laplace transform of discrete truncated power is given by: 

E e-^^-h,[h) = n (2) 
heZ" aex ^ ^ 

That motivates some scholars to study discrete truncated power in another 
way. Szenes and Vergne establish the relation between discrete truncated 
power and Jeffrey-Kirwan residues, and propose a residue method to com- 
pute discrete truncated power Concini and Procesi present the relation 
between discrete truncated power and toric arrangement, and show that dis- 
crete truncated power is a quasi polynomial on some rational lattice [31H] . 

In this paper, we investigate the discrete truncated power by using toric 
arrangement. Our idea comes from two aspects: the relation between mul- 
tivariate spline and hyperplane arrangement |3j, and the Laplace transform 
of the discrete truncated power. Although the second technique was used in 
[T][5] . within the complex field, our argument and result are simple and clear, 
as the discussion is based on the real field. Moreover, we prove that the dis- 
crete truncated power tx is a periodic piecewise polynomial on the integral 
lattice cone spanned by the vectors in X. This is more precise than the re- 
sult of Concini and Procesi in pjHI . Moreover, we propose a toric arrangement 
method to compute the discrete truncated power, which is able to give the 
explicit expression of tx- 

The paper is organized as follows. In Section 2, we briefly review the con- 
clusion given by Concini and Procesi, and introduce some useful definitions 

^ In [3], the authors show that discrete truncated power is a quasi polynomial on the 
rational lattice. However, as the discrete truncated power is defined on the integral 
lattice, their result essentially makes sense merely for the integral lattice. Hence our 
result can be considered as an improvement of theirs. 
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and notations. In Section 3, we propose a new method to show the relation 
between discrete truncated power and toric arrangement, by using factor de- 
compositions in the real field. We can see that our method can be used to 
compute the explicit expression of the discrete truncated power, which is also 
illustrated in two examples. 



2 Toric arrangement and discrete truncated power 

In this section, we introduce some definitions and notations, and previous 
results regarding the relation between toric arrangement and discrete trun- 
cated power. We make the stipulation that in this paper, the set associated 
with the discrete truncated power is multiset. 

Let A be an integral lattice in R'' and 

C[A] := {e<"'^\ a e A} 

be a collection of multivariate exponential functions. Let v be an integral 
vector in and be a translation operator such that 

-Tvfix) = f{x + V) 

for any function / on R*. Then we define 

t[A] = {Ty, veA} 

to be the collection of translation operators associated with A. Denote C, 

to be the operators of the Laplace transform and the inverse Laplace transform 

respectively. Then the elements in C[A] and t[A] have the following relations: 

Let = Xi^X2^ ■ ■ ■ Xg" be a monomial in s variables, where cti, q;2, • • • , cts are 
nonnegative integers, and denote 

pW(a;) := {^aax", a e N", e R} 

a 

to be the space of all the polynomials in s variables. Let D"' = D^^D^^ ■ ■ ■ D"' 
be a difference operators in s variables, and denote 

:= {^a^L*", a e N^ a« e R} 

a 
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to be the space of all difference operators in s variables. Similarly, the elements 
in P^*)(x) and P^'^\D) have the following relations: 

= {-Dr, 
C-\D'') = (-x)". 

Because the Laplace transform of tx'- Eq.([2]), has the form 

1 



p — {a,x) ' 

the following theorem is given. 
Theorem 1 ^ Let 

Sx = C[A]P^'\D) n (1 - e-<"'^>)-\ Tx = r[A]pW(x)tx 

be two collections of functions. Then under the Laplace transform, Tx is 
mapped isomorphically onto Sx- 

This theorem means that for any element Tcp{x)tx{x) in Tx, there exists an 
element e^^'^^p{D) Ylaexi^ ~ ^-ax-^-i such that 

£(r,(p(x)tx(x))) = pix + c)tx{x + c) = e<^'^>p(D) [] (1 - e-<"'">)-^ 

aex 

Contrarily, for any element e^^'^'^p{D) Haexll ~ e^'^^)^^ in Sx-, there exits an 
element Tcp{x)tx{x) in Tx such that 

£-i(e(c,.)pP) -Q (1 _ e-<"'^->)-i) = p{x + c)tx{x + c) = rM^)tx{x)). 

ai^X 



To study the properties of tx on the lattice, Concini and Procesi propose a 
toric reduction method to study the discrete truncated power. The result is 
given in the following theorem. 

Theorem 2"^ Let X = {ai , 02, . . . , ax} be a set of N points on A, and X 
be all the vectors in X with positive rational multiples. Then the function 
Haex can be written as a linear combination of the form WaeA (i_f,-~la.x)-^ 

with constant coefficients, where A is a linearly independent set of elements 
in X . 

By using the preceding theorem, Concini and Procesi give the following theo- 
rem, which describes the general structure of tx- 

[H — 

Theorem 3 Let A be a linearly independent set of elements in X , such 

that the function Yla&x i_e-(a,^) can be written as a linear combination of the 
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form YiaeA (^i^e-(°-'''>)'^o. with constant coefficients. On the closure of the cone 
spanned by A, the partition function tx coincides with a quasi polynomial for 
some rational lattice A/n. 

The idea of Concini and Procesi gives us a new understanding of discrete 
truncated power. Based on this idea, we get our main result in the next section 
(Theorem |6]). 



3 Main results 



We present our results in this section. First of all, we give a useful lemma. 

Let yi, ai, i = 1, . . . , s be 2s functions from to R, such that yi is a 
nonzero function for each i. Then we have the following lemma. 

Lemma 4 Let yo = aiyi + a;2?/2 + ■ ■ ■ + tts^/s be a function and hi be a positive 

s , 

integer, i = 1, . . . , s. Then l/{yo Yi Ui can be written as a linear combination 

s 

of elements of the form f / H Vi j where f is the product of some 1 < 

1=0 

S S 

k < s, ti, nil are integers, and J2 ti = ^ + J2 hi. 

i=0 i=l 



PROOF. Assume that ai is a nonzero function. Then we get the following 
algorithm to decompose l/{yo Vi')'- 



i=l 



am H h asys) «i 



2/0 n y'l' vl n y'l' yl n vT ■ v 

1=1 1=1 '=1 



Then l/(?/o fl vt^^ can be decomposed into a sum of s terms. The numerator 
of each term is In the denominator, there exists y^ whose power decreases 

s 

by one after decomposition, and the sum of powers of all y^ is also 1 + ^hi. 
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By using the algorithm again, we have: 

aj _ aj{aiyi H h asVs 



vlvT^ n i/f' vlv]'-^ n ^ 

i—\ i— 1 



(4) 



voVj' n Vi 1=] vov/ Vk n U. 

■' i=l ■' i=l 



If the denominator of a term in the right hand side of Eq.(jl]) is the form 

s 

of n Ui-, k G {l,...,s}, we stop using the algorithm on this term. Then 

S 

we can express the numerator of this factor in the form of H OiT' t>y using 

1=1 

Eq-Qj where rrii is a nonnegative integer. Because J2i=i hi is a finite number, 

s ^ 

l/(?/o n Vi^) can be decomposed into a linear combination of elements of the 

i=0 

form n aT'/ ft vl', where E = 1 + E 

j=l i = i = i=l 

If some ttj = 0, then we denote cr := {z|aj = 0} to be the collection of 
indices such that a, = 0. As 

1 1 1 



i=0 

We consider the decomposition of l/(yo njG{i,...,s}\(T using the similar 
argument, we can get the result. The proof is finished. □ 



We can see the result of Lemma H] is similar with the reduction of hyperplane 
arrangement. However, our result generalizes from a number to a function. 

Let X C Z'' \ {0} be a set of integer vectors in such that dimX = s. 
Denote X to be the collection of nonnegative integral multiples of the vectors 
in X. Then we can get the following reduction theorem which is different from 
Theorem [2l We shall explain the difference later. 

Theorem 5 The function Yiaex i_e-{a,x} can be written as a combination of 
elements of the form JlaeA (i„e-<°^^))fett ? '^'^^^ ha G N, such that 

Y.ha = #X, 

a&X 

where A is a linearly independent set of s elements in X and c is a linear 
combination of elements in X . 
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PROOF. We use the mathematical induction on the number of elements in 
X. If X is linearly independent, the result is obvious. Otherwise, we assume 
that the conclusion is valid for X with — 1 elements, where N > s + 1. Let 

X = {ai, 02, . . . , ttj^} 

be a set of vectors. As X is not independent, there is an element, say a^, 
which depends on the rest of the system 

X' = {ai,a2, . . ■,a^_^}. 

From the assumption, the function Haex' i_e-ia,x) can be represented as a 
linear combination of elements of the form 



Ha 



e 



(c',x) 



p-{a,x)\ha 



where ha G N, such that Y,a&A' ha = H^X' , d = J2a£X' (here a^can be 
repeated) and A' is a linearly independent set of s elements in X' where 
X' = {na\a G X', n G N}. Then for the function Jlaex i_e-^,^) ' have 



TT t = t V TT - 



where denotes the summation of all the terms with the form OagA' (X ^e-(°'^>)''tt 
Now let us consider the factor 

(1 _ e-K'^)) _ ^-{a,x)Ya ■ 

If the factor ([5]) can be written as the linear combination of the form 



n /I _ p-{a,x)\ha 

where c is a linear combination of elements in X and A is a linearly indepen- 
dent set of s elements in {rial a G X}, n G N}, then the conclusion is valid. 
Without loss of generality, we let A' = {ai, 02, . . . , a^} (then ha = ha^). And 
let hi = ha^ for simplicity. Because 



rs 
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there exists some rrij G N, z G {1, 2, . . . , s, N}, such that 



iecr ie{l,...,s}\CT 



where a C {1, . . . , s}. 



7 



Assume mjT^O, i — 1, . . . , s, N . By reordering the elements in A: 



m 



=k+l 



we have 



k s 
— ]^ g i=l i=k+l 



{niai,x) 



i=l 



s-1 



i=k+l 
k 



i=l i=k+l 1=1 

k-1 

_|_ JJ Q-{n-iai,x)^-^ _ ^-(nkak,x)^ _|_ . . . _|_ (^]^ _ g-("iai.^>^ 



i=l 



s-1 



ris—l 



— g-("iO<>a;) -J-j- ^{niai,x}^_^(neas,x)^^-^ ^-(jas,x)^^-^ _ ^-(ae,x)^ 
i=k+l j=0 



i=l 



_|_ JJ g-K«<>a;)^_g(nfc+iafc+i,a;>-j^ ^ g-(jafc+i,a:)^ (^]^ _ g-(«fc+i,^>^ 



1=1 
fe-1 



j=Q 



"fe-1 



ni — 1 



_|_ JJ g-Kai,a;>^ ^ g-O'sfe-a;)^ _ g-(afc,a;)'^ _|_ . . . -j_ ^ e"^-'"^'^^) (1 — e~^"^'^^). 
j=0 j=0 



i=l 



To give simple argument, we give the following notations: 



ri 



t-l 



0(t = < 



2<t<k + l, 



i=l 
k 



t-l 



i=l i=k+l 
l-l 
j=0 



-{iai,x) 



l<l<k, l^N, 



,{lai,x)Y^^-{iat,x) k + l<l<S, 



1=1 



Vd 



Then we have 



1 - e-^''^'^^ l<d<s, 
1 - e-<'*<^"<^'^> d = AT. 

1=1 



and 

I TT = e^'^''"^— n ^ = fe<^''-)/3^ rr a B-) — 

According to Lemma HI l/(yjv 11^=1 Vi) can be represented as the linear combi- 

s s ^ s i 

nation of elements of the form H oiT /{v^n H yi)i where In+Y, U = hi. 

1=1 i = l «=1 fi = l 

Because each rrii is a nonzero integer in Eq.([n]), 

{Vi-. i = l,...k-l,k + l,...,s,N} 
are linearly independent. We observe that 

(e<'='->/3;vn«,A)(n«r) (7) 
i=i i=i 

is obviously a sum of the form of e^^'^\ where c is a linear combination of 
elements in X. 

If some mi's are equal to in Eq. ([6]). We let 



1=1 



where integer 7^ for z = 1 , . . . ci. Then we consider the decomposition of 

,{c' ,x) 



With the similar discussion, we also can conclude that the hypotheses is right 
for X = {oi, . . . ,a^}. Then the proof is finished. □ 



With Theorem [T] and Theorem |5] in hands, we can get the following theorem, 
which is different from but more precise than Theorem [3l Before proposing the 
theorem, we introduce some definitions. For a matrix X = {ai, . . . , a^} C Z*, 
we denote by 

Ax ■= e N} 

a lattice cone and 

a shift of Ax- Let C{X) be the cone 

C{X) := {^aiGilai G R+}. 
Then the following theorem is established. 



9 



Theorem 6 Let X = {oi, . . . , oat} be a set of integral vectors in 7/ , then tx 
is a sum of some periodic piecewise polynomials on the shifted sublattice cone 
ofAt- 



PROOF. By Theorems [T] and |5l we have 



where A = {ai, . . . , ag} is a basis formed by vectors in X := {na\a E X, n E 
N}, c is a linear combination of elements in X and J^aeA^a = Let 
hi = hai for simplicity. Now, let us compute the inverse Laplace transform of 

Let Af- be a vector who is perpendicular to all the vectors in A except for a^. 
Then we have 



e 



n 



where Da is a difference operator along vector a. Through Eq.(j8]) and Theorem 
[H we can show that the inverse Laplace transform of tt (il'e-(a,x))hi is 

It is clearly that 

^A(a; + XK^fe ~ l)*^/^' + c) = + c) = 1 

k=\ 

on the lattice Aj[(c). Thus the inverse Laplace transform of -pi Q!!e-(a,a:)-)ha 

is a polynomial of degree j^X — s on the lattice Aj^(c). It is obvious that 
Aj[ G C(X) is a sublattice cone of A J. The theorem is proved. □ 



Now let us analysis the difference between Theorem [3] and Theorem [6l The 
main difference is their domain, i.e. the lattice. For discrete truncated power 
tx, the lattice 

Ax = {Sa^ailaj G Z} 

in Theorem [6] is obviously smaller than the lattice Ax /n in Theorem [31 This 
means that tx makes no sense on the point in (Ax/'^)\Ax. The reason comes 
from the different definitions of "X" in two toric reduction theorems: Theorem 
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|5]and Theorem [2l As tx is a function in Z*, and our discussion in Theorem [5] 
is based on R*, while Theorem [2] is discussed in (see |3j). Hence our resuhs: 
Theorems El are the improvements of Theorems HI E] respectively. 

From Theorem El we see that tx is a sum of some periodic piecewise poly- 
nomials on the shifted sublattice cone of Aj^. This property is shown in the 
following two examples. 

Example 7 Let X = {1, 1, 2} be a set of real numbers. Obviously that 

1 



C{t 



X 



(l-e-2^')(l-e" 
As 1 — e""^^ = (1 + — e~^), we have 



1 



(l + e--)2(l 



-x\2 



■x\2' 



xy (1 _e-2x)3(i _e-x)2 (i _ e-2x)3 v ^ 2!(-2)2 1 - e'^^ 



(l-e-2^)(l-e 

5*0 ^/ie inverse Laplace transform of 
(x + 2)(x + 4) 



(1+e- 



(l-e-2^)(l-e-^)2 

= (1 + r-0^ ^-^^T"^' ^2(a: + 4) = (1 + r.,f 



is equal to 

x + 2)(x + 4) 
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t2{x) 



on the lattice cone A^ . Therefore, tx can be expressed in the following simple 
form: 

' 1 X = 

(x + 2)(x + 4) 



tx{,x) 



4 

{x + 2)2 



X e A+(l) 
a; e A+(2) 



Example 8 Let X = {(q), (°), ( 2)} ^ '^^^ of real numbers. Observe that 



C{t 



(1 - e-^)(l - e-!')(l - e-(-^+2J/))^ 



(l + e-2^)(l-e-^') 



-2y -a:+ 



e--)2(l-e-2^)(l-e-(-^+2j/)) 



(1 - e--)2 
1 + 6-^^ 



g-2|/g-a::+2y 



(1 - e-^)2(l - e-(-^-+22')) (1 - e-^)2(l - e-'f) 



(1,0) 



-2 (1 -e-^')(l -e"("^'+22/)) 
where we use the identity 



(1 - e-^')(l - e~^)' 
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1-e 



— X 



(1 _ e-(2|/-(-a;+2j/))^) 

(1 - e-^y) + e-^y{l - e-"+2y) 

(1 + e-^)(l - e-y) - e-22^e-"+2y(i _ e-(--+2j/)). 



Let y4i = {(g), ( 2 )} ffl'^c? = {(g), (^)}, t/icTi t/ic mvcrse Laplace transform 

of (l-e-^)(l-e-i/Kl-e-(-+^^)) ^^^^^ 



= (1 + r(0-l)) 



2x + ?/ + 2 
2 



tAi (x + 1, y) - r(_i,o) (x + l)tA2 (a; + 1, 2/) 



on t/ie lattice cone Aj^ . Therefore, 



2x + y + 2 
2 



tAi(a;,i/) + 



2a; + ?/ + 1 



2 



tAr{x,y- 1) - xtA^{x,y). 
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